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STABILITY  IN  LINEAR  DELAY  EQUATIONS 

by 

Jack  K.  Hale,  Ettore  F.  Infante 
and  Fu-Shiang  Peter  Tsen 


ABSTRACT 

For  linear  autonomous  differential  difference  equations  of  retarded 
or  neutral  type,  necessary  and  sufficient  conditions  are  given  for  the  zero 
solution  to  stable  (hyperbolic)  for  all  values  of  the  delays. 
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1.  Introduction 

This  paper  is  devoted  to  the  study  of  the  effect  of  the  delays  on  the 
asympotic  behavior  of  the  solutions  of  linear  retarded  and  neutral  differential 
difference  equations.  A  special  case  of  the  retarded  equations  considered  is 

(1.1)  *(t)  =  AQx(t)  +  Akx(t-rk) 

where  X  €  Fn,  each  A^  is  an  n  x  n  matrix  and  each  r^  >  0,  k  *  1,2,...,N. 
It  is  known  that  the  asymptotic  behavior  of  the  solutions  is  determined  from 
the  solutions  of  the  characteristic  equation, 

(1.2)  f(X,r,A)  dif  det [XI  -An-  f  A.e  _Xr*]  -  0, 

u  k*l  K 

Let  o(r,A)  be  the  supremum  of  the  real  parts  of  the  X  satisfying  (1.2). 

It  is  well-known  that  a(r,A)  <  0  implies  the  zero  solution  of  (1.1)  is  uniformly 
asymptotically  stable  (see,  for  example,  [4]). 

Because  the  supremum  o(r,A)  is  attained  at  some  specific  value  of  X 
satisfying  (1.2)  and  the  function  f(X,r,A)  is  continuous  in  r,A,  it  follows 
that  a(r,A)  is  continuous  in  r,A.  Therefore,  the  property  of  being  asympto¬ 
tically  stable  at  some  point  r°,A°  is  preserved  under  small  perturbations 
in  r,A  from  r°,A°. 

Our  primary  objective  is  to  give  conditions  on  the  coefficients  A  in 
(1.1)  which  will  ensure  that  Eq.  (1.1)  is  asymptotically  stable  for  all  delays 
r  ■  (r^,  ...  rN)  with  rfc  >  0,  k  ■  1,2,  ...  ,  N.  that  is,  we  want  to  charac¬ 
terize  those  values  of  A  such  that  o(r,A)  <0  for  all  rk>0,  k  *  1,2,  ..,  N 
Some  aspects  of  this  problem  have  been  previously  discussed  by  Zivotovski  [7] , 
Dstko  [3],  Repin  [5],  Silkowskii  [6],  Cooke  and  Ferreira  [2]. 
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If  we  call  the  set  of  such  A  the  stable  cone  S  for  (1.1) ,  then  one 
of  our  results  states  that  ACS  if  and  only  if 

(i)  det  lU  A.  t  0 

k-0  K 

(ii)  det [iy  -  K,  -  f  A.**1  ....  sN]  i  0  fbr  all 

°  k-1  K  1  ff 

y  €  R,  y  *  0.  sj  €  C,  |Sj|  -  1,  j  -  1.  2,  ...  ,  N. 

rll 

These  conditions  also  iaply  that  the  spectrua  of  V  A.  and  the  spectrua 

k-0 

of  Aq  lie  in  the  left  half  plane. 

If  the  equation  (1.1)  is  a  scalar  equation,  then  the  above  conditions 
for  A  to  be  in  S  simplify  to  .  A^  <  0,  J**  jA^j  |Aq|.  letter 

result  was  obtained  by  Zivotovski  [7]. 

In  the  applications,  it  is  not  always  true  that  the  delays  vary 
independently  of  each  other.  For  exaaple,  with  three  delays,  ri»r2*r3»  on( 
nay  have  Xj  ■  Sj,  r2  »  s2,  «  Sj  ♦  s2  for  some  positive  nuabers  s|*s2*  In 

this  case,  the  stable  cone  can  be  larger  than  the  one  obtained  before.  Me  also 
give  a  characterization  of  the  stable  cone  in  this  case. 

Finally,  the  results  are  extended  to  the  such  no  re  complicated  case  of  a 
neutral  differential  differente  equation 

•**■<*>  -  -  Ao  * 

The  basic  difficulty  here  arises  from  the  fact  that  the  asymptotic  behavior  of 
the  solutions  of  the  difference  equation 

x(t)  -  I*  i.  x(t  -  r.)  -  0 

k-1  *  * 

depends  in  a  very  conplicated  way  upon  the  delays  rfc.  The  results  far  tills 
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case  are  contained  in  Section  5  and  rely  heavily  upon  Avellar  and  Hale  [1]. 

2.  General  results  for  retarded  equations. 

Suppose  •  [0,*),  R  ■  (-",»),  r  ■  (tj,  €  QR  )  , 

Yk  -  (Ykl.  »  YkM).  Ykj  >  0  integers,  Yk  t  0,  Yk*r  *  y^Ty 
k»l,2,  . .  ,  N,  j»l,2,  ..  ,  M,  and  consider  the  retarded  differential  difference 

equation 

N 

(2.1)  x(t)  -  AQx(t)  ♦  £  Akx(t  -  Yk*  r) 

where  x  €  Rn  and  each  A^,  k«0,l,  ...»  N  is  an  n  x  n  real  constant  aatrix. 
The  characteristic  function  for  Eq.  (2.1)  is 

_«  ~xYk-r 

(2.2)  f(A,r,A)  -  detlAI-AQ-  T  Ake  1 

k«l 

2 

where  A  -  (Ap.Aj,  •••»  \)  €  lRn  . 

Definition  2.1.  System  (2.1)  is  said  to  hyperbolic  at  (r,A)  if  f(\,r,A)  *  0 
i up lies  Re  A  i  0.  System  (2.1)  is  said  to  be  asymptotically  stable  at  (r,A) 
if  f(A,r,A)  *  0  implies  Re  X  <  0. 

The  delays  in  Eq.  (2.1)  are  the  constants  Yk'r»  h*l,2,  ...»  N.  They  are 
not  independent  and  are  determined  by  the  vector  r  »  (r^ , . . .  ,rM)€(JR+)M.  For 
example,  if  M  -  2,  N  -  3,  r  •  (r^),  y^  (1,0),  y2*  (<M).  Ys  »  (1.1). 
the  delays  are  r1,r2,rl  +  x^. 

Our  objective  is  to  determine  conditions  on  the  coefficients  A  in  (1.1) 
to  ensure  stability  (or  hyperbolicity)  for  all  values  of  r€(]R+)*®.  This  means, 
in  particular,  that,  for  a  given  r°,  we  must  have  stability  (or  hyperbolicity) 
for  all  ar°  with  a  >  0.  By  letting  t  •*>  at  in  (1.1),  this  means  that,  if 
A0  ensures  stability  (or  hyperbolicity)  for  all  r€  (R+)M,  then  a  A®  also 
ensures  this  for  every  o  >  0;  that  is,  the  set  of  such  A’s  is  a  cone.  No 
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foraalize  these  ideas  in  the  following 

+  M 

Definition  2.2.  For  a  given  r£(R)  ,  the  ray  Yr  through  r  is  the  set 
{a  r  a  ^  0).  For  a  given  r°  €(*+  )**,  the  hyperbolic  cone  at  r°, 

designated  by  H^o,  is  defined  by 

2 

H  »  U6Sn  .  Eq.(2.1)  is  hyperbolic  at  (r,A) 
r° 

for  every  r  €  y  Q) 

T 

The  hyperbolic  cone  H  is  defined  by 

H  =  n(Hr:r€(R+^) 

For  a  given  r°€(!R+/*,  the  asymptotically  stable  cone  at  r°,  designated  by 
S  D,  is  defined  by 

2 

S^o  *  {A€Rn  ^N+1^:Eq.  (2.1)  is  asymptotically 
stable  at  (r,A)  for  every  r  €  Yr°} 

The  asymptotically  stable  cone  S  is  defined  by 

s  »  nr(Sr:r  €(K+)M) 

In  the  following,  the  notation  ReX(A)  for  a  matrix  A  designates  the  set 

consisting  of the  real  parts  of  the  eigenvalues  of  A.  As  a  preliminary  for 

the  classification  of  Hr,Sr,  we  have  the  following  elementary  result,  a  form 

of  which  was  proved  by  Datko  [3]. 

Theorem  2.3.  A  €  H  [or  SI  if  and  wily  if 

*  •  » 

N  N 

(i)  Re  X  (I  A. )  j*  0  [or  Re  X(  \  A. )  <  0] 

k-0  *  ■“*  k-0  * 

(ii)  f(iy,ar.  A)  i  0  for  all  y  €  R,  y  +  0,  o  >  0. 
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PTgof:  Let  so  ■  max  {Re  X£0:f(X,  ar„A)  »  0>,  ua  ■  min  {Re  X  >  0:f(X,ar,A)  •  0} 
with  s  ■  - u  if  the  corresponding  set  is  empty.  The  numbers  s  ,u 

are  continuous  in  a.  Condition  (i)  implies  sQ  <  0,  Ug  >  0.  If  AftHr, 

then  there  is  an  ot>  0  such  that  either  *  0  or  u  »  o.  This  con- 

0  ao  ao 

tradicts  (ii).  Thus,  (i),(ii)  imply  A€Hr.  The  converse  is  obvious. 

For  the  more  difficult  equations  with  distributed  delays,  Cooke  and 
Ferreira  [2]  have  obtained  nontrivial  results  in  the  spirit  of  Theorem  2.3. 

Theorem  2.4.  A  €  H  if  and  only  if 
N 

(H.)  det  l  A.  f  0 
1  k*0  K 

Jj  y  y 

(H2)  det  [iy  -  AQ  -  l  A|c5ikl*--sMklB^  *  0  for  ali  y£*»  *  *  0 

s.€  C,  | s . {  =  1,  j  *  1,2,  ...  ,  M. 

3  J  N 

ACS  if  and  only  if  (Hj)  and  ReX(J  A^)  <  0. 

Proof.  Suppose  A€H.  Then  (Hj)  is  satisfied.  If  (H2)  is  not  satisfied 
at  (y.s),  choose  y  /  00^,  so  that  -yS^  >  0  and  s^  «  exp(i0k)  for  all  k.  With 

a  -0k/y,  we  have  f(iy,r,A)  *  0  which  contradicts  the  fact  that  A  €  H. 

Conversely,  suppose  (H^) ,  (H2)  are  satisfied  and  A  f.  H.  Then  there 

is  an  o  >  0,  y  €  R,  y  M,  r°  €  (»+)M  such  that  f(iy,  ar°.  A)  »  0.  Since 

this  contradicts  (H_) ,  we  have  proved  the  first  part  of  the  theorem. 

2  N 

The  condition  ReX(J  A. )  <  0  is  equivalent  to  saying  that  Bq.(2.1) 

k-0  * 

is  asymptotically  stable  for  r  -  0.  Thus,  the  last  statement  in  the  theorem 
is  true. 

To  obtain  other  characterizations  of  H,S,  we  need  the  following  1mm- 
We  are  grateful  to  John  Mai let- Par et  for  assistance  in  the  statements  and 
proofs  of  the  next  two  results. 
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La— a  2.S.  .  If 

N  Ykl  Y;:J, 

(2.3)  P(X,Sj,  ...  ,  s)p  *  dot  j*k*l  J* 

then  the  hypothesis  (H2)  iap lies 

(2.4)  P(iy»slt  ...  ,s^  p  0  for  y£  R,  |s^|<  1.  j  -  1,2,  ...  ,  M. 
Proof:  Fix  s®,  |s?|  «  1,  j  -  1,2,  ...  ,  M  end  consider  the  function 

Q(X,a)  *8*  P(X,asJ,  ....  osjj). 


Designate  the  zeros  of  this  equation  by  a(X) .  Then  a(X)  is  —ronorphic, 
defined  on  sone  Reinann  surface  over  the  A-plane  and  lia  jx|-^c»  |a(X)|  «  -. 

Consider  the  curve  in  the  a-plane  defined  by  a(iy) ,  y  €  ».  Hypothesis  (H2) 
implies  that  |a(iy)|  >  1  if  y  p  0.  Thus,  |a(iy)|  >  1  for  all  y.  Thus. 
Q(iy.a)  -  P(iy,asJ,  ...  ,  a$  P  0  for  |a|  <  1,  y  €R  and  all  s® •  1  s® |  «  1, 


j  *  1,2,  ...  ,  M.  This  proves  the  le— a. 

By  taking  each  s^  »  0,  j  ■  1,2,  ...  ,  M, in  Le— a  2.5,  we  obtain 

Corollary  2.6.  Hypothesis  (Hj)  iaplies  Re  X  (Aq)  p  0. 

N 

Corollary  2.7.  If  A  €  H,  then  Re  X  (J  A*.)  p  0,  Re  X  (AfP  p  0  and  the 
N  *"° 

—trices  l  A^  and  A,,  have  the  sa—  number  of  eigenvalues  with  positive 
k*0  * 

and  negative  real  parts. 

If  ACS,  then  Re  X  0U  <  0,  Re  X  (J  AJ  <  0. 

-  u  k«0 

Proof;  Let  P(X,s)  be  defined  by  Relation  (2.5).  Let  Q(X,s)  ■  P(X,  s,  ...  ,  s). 
Ft—  Le— a  2.5  and  Relation  (2.4),  Q(iy,a)  p  0  for  0<_s  <  l.y  €*.  Hypothesis 
(H1)  iaplies  Q(iy,  1)  p  0  for  all  y  €*.  Thus,  Q(- ^)  has  no  roots  on  the 
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imaginary  axis  for  any  s  6  [0,1].  Since  Q(X,0)  *  det  [XI  -  Afl],  Q(X,1)  * 

N  rN 

det  [XI  -  l  A^] ,  it  follows  that  Re  X(AQ)  +  0,  ReXQ,  A^  i  0.  Since  the 


k-0  "■  w  k=0 

zeros  of  Q  (X^)  »  0  are  continuous  in  s,  the  result  follows  immediately 


The  verification  of  Hypothesis  (H^)  is  extremely  difficult.  For  the 


case  of  an  n^  order  scalar  equation  and  independent  delays,  this  hypothesis 
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i 

f 

l 


l 


■  *>■ 


Proof:  Transform  the  equation  to  an  equivalent  system 

N 

x  =  AQx(t)  +  J  iArx(t  ‘ 

where  x  =  (y,y^\ . . .  ,y^n_1^)  .  The  characteristic  function  is  given  by 

N 

f(X,r,a)  =  pQ(X,a)  +  £  ^(X.aje-^ 

Hypothesis  (H1)  is  equivalent  to  f(0,r,a)  t  0  which  is  (2.6).  Hypothesis  (H2) 
is  equivalent  to 

pQ(iy,a)  +  l  pk(iy,a)  sk  /  0  V  y  €  1R,  y  ^  0,  | skl  =  1, 

and  this  is  equivalent  to  (2.7).  The  last  statement  follows  from  Corollary  2.7 
since  the  characteristic  function  for  AQ  is  pQ(X,a).  This  proves  the  theorem. 

We  also  can  generalize  Zivotovskii's  result  to  the  case  where  the  delays 
are  dependent,  but  it  cannot  be  stated  in  such  a  simple  fashion.  The  proof  is 
the  same  as  before. 

Theorem  2.9.  Consider  again  Eq.  (2.5)  with  uk  =  Yj*r,  r  €  (R  )  ,  Yj  *  (Yjj* • • • *YjM) » 
Y ..  nonnegative  integers,  y^  /  0.  With  Pq(X»*)»  Pk(X»a)  defined  as  in  Theorem 

JK  J 

2.8,  the  vector  a  €  H  if  and  only  if  Relation  (2.6)  and 

N  Yk i  Yku  . 

(2.9)  pQ(iy,a)  +  J  ^Pk(iy,a)  Sj  *  0  Vy€F,  y  +  0,  |s^|=  1, 

j  *  1 1 2 1 « •  • 

are  satisfied. 

The  vector  a£S  if  and  only  if  relations  (2.6)  ,(2.9)  ,(2.8)  are  satisfied. 

For  the  case  of  one  delay  in  Eq.  (2.1)  one  can  obtain  an  equivalent 
formulation  of  the  cones  H,S  following  an  idea  of  Repin  [5]  via  the  following 
lemma. 
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Leana  2.10.  If  Ag*Al»  are  real  n  x  n  matrices ,  then  the  statement 
(2.10)  For  every  y  €  1R,  y  t  0,  the  solutions  of  the  equation 


det 


AAq+Aj 

-Xyl 

Xyl 

XAq+A 

satisfy  | X|  <  1 
is  equivalent  to  the  statement 


(2.11)  For  every  y  £»,  |a|  <  1,  det[iy  -  AQ  -  aA^  t  0 

Proof:  If  X  /  0,  y  i  0,  then  the  equation  in  (2.10)  is  equivalent  to 

'  0  -Xwl  .  (XAq.A,)2 

det  -  0 
Xyl  XAq+Aj 


which  is  equivalent  to 

det  (X2y2I  +  (XAq+Aj)2]  -  0 

which  is  equivalent  to 


det  [-X(iy)  +  (aAq+Aj)]  *  0 

for  y  €  R,  y  i  0.  If  [X|  <  1  when  this  is  satisfied,  then  this  is  equivalent 
to  saying  that 


det  [iy  -  AQ  -  aAj]  -  0,  y  €  R,  y  +  0 

implies  |a|  >  1.  Thus,  the  solutions  of  this  equation  for  all  y  CK  satisfy 

|o|  1.  Thus,  we  obtain  (2.10)  is  equivalent  to  (2.11)  and  the  lemma  is  proved. 
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Corollary  2.11.  For  the  n-diaensional  system 

*( t)  **  AgX(t)  ♦  AjX(t-r) 
A  ■  (Ag,Aj)  6  H  if  and  only  if  (2.10)  and 


(2.12) 


det(AQ+A1)  +  0 


are  satisfied.  A  €  S  if  and  only  if  these  conditions  and  Re  X (A^)  <  0 . 
Proof:  This  is  an  immediate  consequence  of  Lemmas  2.10,2.5  and  Theorem  2.4. 

3.  First  order  scalar  equations.  For  first  order  scalar  equations, 

N 

(3.1)  *(t)  ■  aQx(t)  +  \  •kx(t-YkT) 

where  a  £F,  j  ■  0,1,..., N,  the  characterization  of  the  hyperbolic  and  asymp¬ 
totically  stable  cones  can  be  specified  in  terms  of  properties  of  the  solutions 
of  the  difference  equation 


N 

(3.2)  aQ +  X  *  °» 

It  is  the  purpose  of  this  section  to  obtain  such  a  characterization. 

Let  a  »  (a0,a1,...tak) , 

N 


(3.3) 


o(0,a)  *  a Q  *  l  a.  cos  Yk’0 

u  k*i  *  * 

N  y 

8(0, a)  •  l  a.  sin  Yk*0  ,  0  €*" 

k-1  *  * 


The  characteristic  equation  for  the  difference  equation  (3.2)  is 

.  *  -Xy,,.* 

g(X,a,r)  d8*  a-.  ♦  J  a.e  k 
w  k«l  * 


(3.4) 
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The  functions  a(0,a) ,  8 (6, a)  are  related  to  the  function  g(X,a,r)  by 
the  relation 

(3.5)  g(iy,a,r)  ■  a(-yr,a)  -  i8(yr,a) 

The  main  result  of  this  section  is  the  following 

Theorem  3.1.  For  the  scalar  equation  (3.1),  a€H  if  and  only  if 
N 

(H„)  l  a.  f  0 
k*0  K 

(H4)  For  each  8  either  a(6,a)  i  0  or  simultaneously, 

a(8,a)  ■  0,  8 (8, a)  *  0.  If  (H^)  is  satisfied .then  aQ  i  0. 

The  vector  a  €  S  if  and  only  if  (Hj) , (H^)  and  ap  <  0 . 

Proof:  For  n  *  1,  the  condition  (Hj)  in  Theorem  2.4  is  equivalent  to  (Hj) 
and  (H2)  is  the  same  as 

(a(8,a),  y  -  B(8,a))  f  0  for  all  y€  R,  y  i  0,  8  €RM 

which  is  equivalent  to  (N^).  Corollary  2.7  implies  aQ  f  0.  The  statement 
about  S  is  also  a  consequence  of  Corollary  2.7.  This  completes  the  proof 
of  the  theorem. 

An  immediate  consequence  of  Theorem  3.1  is  the  following  result  of 
Silkowski  [8]. 

Corollary  3.2  Suppose  the  components  of  r  are  rational  and  define  the 
functions 

N 

Y  (y)  »  «o  ♦  I  ak  cos  yk  •  r  y 
N 

$  (y)  •  I  «L  *in  Yk  •  r  y 

k«l  *  * 

Then  a  €  Hr  if  and  only  if 
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(i) 


“k*0 


(ii)  For  each  y  €  F,  y  +  0,  either  y(y)  t  0  or  y(y)  *  0,  6 (y)  »  0. 
Furthermore,  a  6  Sr  if  and  only  if  (i),(ii)  and  aQ  <  0.  Finally,  condition 
(ii)  is  equivalent  to 

Re  X  <  0 

Proof;  Since  the  components  of 
has  only  one  independent  delay. 

It  is  interesting  to  state  Corollary  3.2  in  terms  of  properties  of  zeros 
of  polynomials.  If  the  components  of  r  are  rational, we  can  write 


**  -*Vr 

if  an  +  i  a.  e  «  0. 

—  u  k«l  K 


r  are  rational,  we  may  assume  the  equation 
The  result  is  then  a  special  case  of  Theorem  3.1 


N 

^  1  ^  c°s  Yk  ■  r  y  *  h(cos  y) 


(3.8) 


N 

l  a.  sin  y.  •  r  y  *  (sin  y)  g  (cos  y) 
k*l 

where  h,g  are  polynomials. 

Corollary  3.3.  With  h,g  as  in  relation  (3.8) .the  statement 

(i)  h(i)  i  0  and,  for  every  n  €  [0,1)  for  which  h(n)  *  0, 


it  follows  that  g(n)  =  0 
is  equivalent  to  the  statement 

N  -\y.  *r 

(ii)  h(l)  +  0,  a.  +  0,  Re  X  <  0  if  an  ♦  l  a.e  *  -  0 

u  -  0  k»l  * 


Proof;  This  is  a  restatement  of  Corollary  3.1. 
Corollary  3.4.  For  the  scalar  equation 

N 

*(t)  -  aQx(t)  ♦  l  akx(t-rfc) 
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N  N 

we  have  a  €  H  if  and  only  if  I  a».  i  0,  I  |aj  <  |aj . 

k«0  K  k-1 

Tha  vector  a  €  S  if  and  only  if  a  €  H  and  a^  <  0. 

N 

Proof:  Theorea  3.1  iaplies  a  €  H  if  and  only  if  Bq  t  0,  J  a.  4  0  and 

k*Q 

(H4).  Hypothesis  (H4)  is  equivalent  to  either 


N 


a(Q.a) 


viw\*«V0 


or,  siaultaneously. 


o(0, a)  ®  0 
N 

S(6, a)  «  l  a.  sin  0.  »  0. 
k-1 
N 

Tha  latter  relation  iaplies  £  |a.|  <  |a_|.  Conversely,  if 

k-1 

N  rN 

l  | a.  |  <  |aJ,  then  o(0.a)  «  0  for  some  0  iaplies  \  |a.|  »  jaJ 

k-1  K  ~  u  k-1 

and  each  component  0^  of  0  is  0  or  it.  But  this  iaplies  0(0, a)  *  0. 

The  last  assertion  about  S  is  also  a  consequence  of  Theorem  3.1. 

Corollary  3.5.  For  the  scalar  equation 


x(t)  *  aQx(t)  ♦  a1x(t-r1)  ♦  a2x(t-r2)  ♦  a^t-rj-r^ 
we  have  a  €  H  if  and  only  if 


(i)  a0  ♦  ♦  a2  +  a3  i  0 


(ii)  1  ♦ 


0  < 
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(iii)  either  ^  ^  *  a3  or  '  *3* 

The  vector  a  €  S  if  and  only  if  a  €  H  and  aQ  <  0. 

Proof:  Theorem  3.1  iaplies  that  a  €  H  if  and  only  if 

aQ  i  0,  (i)  and  (H^)  is  satisfied;  that  is,  if  ■  a^  /  a Q,  and 

a(8,a)  «  aQ  {1  ♦  bj  cos  0j  ♦  b2  cos  02  ♦  bj  cos  (Oj+Oj)], 

6(9, a)  =»  aQ  [bj  sin  0j  ♦  b2  sin  02  ♦  bj  sin  (0j  ♦  02)]» 

then,  if  there  is  a  0  such  that  a(0,a)  *  0,  then  6(0, a)  »  0. 

But,  these  0  are  precisely  the  ones  for  which  a(0,a)  ♦  i8(0,a)  -  0; 
that  is. 


i0. 


i0„ 


ie. 


1  ♦ 


bjC 


-e 


fb2  +  b3® 


To  have  a  solution  of  this  equation,  one  Bust  have  0j,02  take  on  tb*  values 
0  or  it.  Checking  separately  the  cases  ®i  “  ®i  “  11  *  one  observes  that, 
for  any  b^,  one  has  a  solution  of  this  equation  if  and  only  if  either 


(3.9) 


1  ♦  b. 


b2*b3 


or 


(3.10)  1  -  bj  ■  b2  -  bj 

If  we  choose  a  value  of  b  ■  (b^.b^bg)  which  does  not  satisfy  either  of 
these  inequalities  and,  if  (H^)  is  satisfied,  then  we  aust  have  a(0,a)  i  0 

for  all  0  »  (0j,  02)  €R2.  The  relation  o(0,a)  f  0  for  all  0  €R2  is 

equivalent  to 

i0,  i09  i0, 

Re  U  ♦  bje  4)  ♦  Re  e  *fl>2  ♦  e  l]  i  0 


for  all  6  €»2.  This  implies  |bj|  <  1.  If  jbjl  <  1,  then  it  is  easily 
observed  that 

if  0  <  bj  <  1,  then  1  -  b^  >  |bj  -  bjj, 

(3.11) 

if  0  <-  bj<  1,  then  1  ♦  b^  >  |b2  ♦  bj| . 

Relations  (3.9),  (3.10)  and  (3.11)  are  precisley  the  relations  in  (ii). 
Thus,  (H4)  implies  (ii) .  The  converse  is  a  straightforward  reversal  of  the 
argument.  This  proves  the  corollary. 

If  there  are  three  independent  delays  in  the  equation  in  Corollary  3.5, 
the  condition  (ii)  would  be  replaced  by  | |  +  |a2|  +  |a3|  <_  |aQ| 
which  is  a  more  restrictive  condition  on  the  coefficients  than  the  one  for 
only  two  independent  delays. 

4.  Some  examples.  In  this  section,  we  give  some  examples  illustrating  the 
application  of  the  results  of  Section  2  to  equations  of  order  >  2.  These 
examples  will  also  show  that  the  results  in  Section  3  do  not  generalize  to 
sy st eats;  that  is,  one  cannot  reduce  the  discussion  of  the  hyperbolic  and 
stable  cones  to  the  discussion  of  properties  of  difference  equations. 

Example  4.1.  Consider  the  system 

(4.1)  *(t)  -  B[x(t)  -  yx(t  -  r)J 

where  B  is  a  2x2  matrix  with  ReX(B)  <0,  u  is  a 
scalar,  |y|<  1.  be  want  to  determine  conditions  on  B,y  so  that  the 


itrices  (B,  -y  B)  €  S,  the  asymptotically  stable  cone. 
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By  us lag  the  Jordan  normal  form  for  I,  one  sees  that  Theorem  2.4  asserts 
that  (B,-  u  B)  €  S  if  and  only  if 

iy  -  X(B)  a  ♦  tie16)  #0  for  all  y  €  *,  y  f  0,  0  €  *, 


for  every  eigenvalue  X(B)  of  B.  If  X(B)  is  real,  this  relation  is 
always  satisfied.  However,  if  X(B)  is  complex,  this  may  not  be  true.  In 
fact. 


1  +  ye 


i9 


2  1/2 

Cl  ♦  2y  cos  0  ♦  y  )  exp  i  t(y,6) 


where 


tan"1  s(y,0)  «(ysinB)Al  ♦  y  cos  6) 

«*»<*  0  <  U(y.B)  I  <  C(u»0)  «  0.  If  CQ(y)  »  naxejc(y,e) | ,  then 

CQ(y)  >  w/2  if  0  <  y  <  1,  C^(y)  <  ir/2  if  -1  <  y  ^  0. 

We  can  now  assert  that  (B,-yB)  €  S  if  and  only  if 


-y-  <  arg  X(B)  ♦  C0(V»)  <  it 

for  all  eigenvalues  X(B)  of  B.  This  clearly  puts  a  restriction  on  the 
eigenvalues  of  B  and  y  with  the  restrictions  being  more  severe  for 
y  >  0  than  for  u  <  0. 

This  example  shows  that  the  results  in  Section  3  cannot  be  generalised 
to  systems.  In  fact,  the  zero  solution  of  the  difference  equation 

•y(t)  -  yiy(t-r)  -  0 

is  asymptotically  stable  for  every  y,  |y|  <  1,  and  (B,-yS)  may  not  belong 


to  S  for  every  B,y . 

Example  4.2.  Consider  the  eqviation 
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(4.2)  X(t)  +  aQ*(t)  +  aj*(t-r)  ♦  a2x(t)  ♦  a3x(t-r)  -  0 
along  with  the  characteristic  equation 

X2  ♦  X  (a^  ♦  aje"*r)  +  a2+a3e”*r  m  f| 

If  a  *  (a0,a1,a2,aJ) ,  then  Theorem  2.8  implies  that  a  €  H,  the  hypexbolic 
cone,  if  and  only  if 


(4.3) 

(4.4) 


a2  +  a3  *  0 


| P(iy) |  >| Q(iy)|  for  all  y  t  0,  y  €F 
P(X)  *  X2  ♦  aQX  +  a2 


Q(X)  =  ajX  +  a3 

The  condition  (4.4)  is  equivalent  to 

f(X2)  d*f  (y2-«2)2  ♦  (a2-a2)y2-a2  >  0. 

for  all  y  €  *,  or,  equivalently,  f(r)  >  0  for  r  >  0.  It  is  easy  to  show  the 
quadratic  function  f(r)  >  0  for  r  >  0  if  and  only  if 


(4.5) 


a2  -  a2  -  2a2  >  0  implies  |a2|  >  |a3| 
a0  '  al  ‘  2a2  <  0  i*^lla*  a2  >  a3  *  ^“O**!*2*^ 


Thus,  a  €  H  if  and  only  if  (4.3),(4.S)  are  satisfied. 


i 
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Theorem  2.8  also  implies  that  a  6  S,  the  asymptotically  stable  cone, 
if  and  only  if  (4. 3), (4. 5)  and  Re  X  <0  if  X^  ♦  a^  X  ♦  a2  *  that  is, 

(4.6)  aQ  >  0,  a2  >  0. 

If  Eq.  (4.2)  is  transformed  to  a  system  of  order  two 

*(t)  -  AQx(t)  ♦  AjX(t-r) 

then  the  zero  solution  of  the  difference  equation 

AQy(t)  Ajy(t-r)  ■  0 


is  asymptotically  stable  if  |a3|  <  |a2|  which  does  not  imply  anything 
about  H. 

S.  Neutral  equations.  In  this  section,  we  generalize  the  results  of 
Section  2  and  3  to  neutral  differential  difference  equations 


(5.1) 


4r [x(t)  - 


N 

*  Aq  x(t)  +  J  x(t-Yk  r) 


where  x  €  Fn,  each  k»l,2,...,N,  is  an  n  xn  constant  matrix 

and  the  yk»r  are  the  same  as  before.  The  characteristic  function  for 
Eq.  (5.1)  is 


(5.2) 


N 

g(X,r,A,B)  -  det  X(I-E  \ 

k*»l  * 


e 


>  r  N  ~Xy 

)-*o-I  v 

k«l  * 


where  A  *  (Ag,A^,. . . ,A^) ,  B  *  (Bj , . > . ,B^) • 

Definition  5.1.  System  (5.1)  is  said  to  be  hyperbolic  at  (r,A,B)  if  there 
is  a  6  >  0  such  that  (ReX:g(X,r,A,B)  ■  0)  n  [-6  ,5]  ■  0. 

System  (5.1)  is  said  to  be  (uniformly)  asymptotically  stable  at  (r,A,B)  if  it 
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it  hyperbolic  at  (r,A,B)  and  {ReA:g(A,r,A,B)  ■  0}n[-S,»)  «  0. 

Definition  S.2.  For  a  given  r®€  CR+)M,  the  hyperbolic  cone  at  r\  designated 

by  H  Q,  is  defined  as 
ru 

2  2 

H  o  *  {(A,B)  €Rn  ^xln  N  Eq.(5.1)  is  hyperbolic  at  (r,A,B)  for  every 
r 

r  ■  ar°,  a  0} 

The  hyperbolic  cone  H  is  defined  by 

H  -  n  (Hr:r€  (**)M> 

For  a  given  r°  €  (R+)M,  the  asymptotically  stable  cone  at  r°,  designated  by 

S  is  the  set  of  (A,B)  €  H  Q  such  that  (r,A,B)  is  asyaptotically  stable 

r  r 

for  every  r  *  ar°,  a  0.  The  asyaptotically  stable  cone  S  is  defined  by 

S  »  nr{Sr:r€0R+)M} 

The  set  Hr  is  not  really  a  cone  in  (A,B)  space.  In  fact,  if  (A,B)  €  Hr, 
then  (aA,B)  C  for  every  a  0,  but  (aA,  oB)  will  generally  not  be.  The 

reason  for  this  is  that,  if  r  **■  ar,  t  ^  at,  then  the  new  equation  has  co¬ 

efficients  (oA.B).  In  spite  of  this  fact,  we  retain  the  tern  cone  for  Hr, 
but  it  should  be  reaenbered  that  the  property  of  being  a  cone  holds  only  in  the 
A  variable. 

Our  objective  is  to  give  a  classification  of  the  hyperbolic  and  asyapto¬ 
tically  stable  cones.  This  problea  is  auch  sore  difficult  than  the  corresponding 
one  for  the  retarded  equation  in  Section  2  because  the  set  of  real  parts  of  the 
zeros  of  the  characteristic  function 

(5.3)  e(A,r,B)  -  det|l  -  J* 

I  k-1 


4 


-20- 


of  the  difference  equation 

H 

(5.4)  y(t)  -  l  *ky(t-yk-r)  -  0 

does  not  depend  continuously  on  r. 

For  the  difference  equation  (5.4), we  need  the  definitions  an&lgous  to 
Definitions  5.1  and  5.2  for  Eq.  (5.1).  '  ■ 

Definition  5.3.  System  (5.4)  is  said  to  be  hyperbolic  at  (r,B)  if  there  is  a 
6  >  0  such  that  (ReX:e(X,r,B)  *0}  n  [-6  ,5]  -  0.  System  (5.4)  is  said  to  be 
(uniformly)  asymptotically  stable  at  (r,B)  if  it  is  hyperbolic  and 
(ReX:e(X,r,B)  «  0)  0  (-6,*)  *  0. 

If  Equation  (5.4)  is  hyperbolic  (asymptotically  stable)  at  (r,B) ,  then 
it  is  hyperbolic  (asymptotically  stable)  at  (ar,B)  for  every  a  >  0.  If  we 
assume 

N 

(5.5)  det(I  -  l  BJ  i  0, 

k«0 

that  it  is  also  hyperbolic  at  o  ■  0  because  (ReX:e(X,0,B)  •  0}  is  empty. 
Thus,  with  (5.5),  if  Bq.  (5.4)  is  hyperbolic  at  (r°,B)  it  is  hyperbolic  for 
every  (ctr,B),  ot  >  0.  This  means  there  is  no  reason  to  use  the  concept  of 
hyperbolic  cone  at  r°  for  Eq.  (5.4).  It  becomes  only  necessary  to  discuss 
whether  or  not  hyperbolic  is  or  is  not  preserved  under  variations  in  r. 

Throughout  this  section,  we  assume  (5.5)  is  always  satisfied. 

Definition  5.4.  Bq.  (5,4)  is  said  to  be  hyperbolic  locally  at  (r° ,B)  if  there 
is  a  mel^horhood  U(r°)  of  r°  such  that  Bq.  (5.4)  is  hyperbolic  at  (r,B) 
for  every  r  C  U(r°).  Eq.  (5.4)  is  said  to  he  hyperbolic  globally  at  ,rB  -  if 
it  is  hyperbolic  at  (r,B)  for  every  r  €  (R+)H.  The  hyperbolic  set  for 
Bq.  (5.4)  is  the  set  (B  C  Hft2|,:Bq.C5.4)  is  hyperbolic  et  every  (r,B)  ,r  € 


til1 


Similar  definitions  are  made  for  asymptotically  stable  locally  at  (r  ,B) , 
asymptotically  stable  globally  at  B  and  the  asymptotically  stable  set  for 
Eq.  (S.4). 

In  Definition  5.4,  no  mention  is  made  of  the  variation  of  the  concept 
of  hyperbolic  with  respect  to  variations  in  the  coefficients  B.  The  reason 
for  this  is  that  if  Eq.  (5.4)  is  hyperbolic  at  (r®,B®),  then  there  is  a  neigh¬ 
borhood  V(B°)  such  that  Eq.  (5.4)  is  hyperbolic  at  (r°,B)  for  every 
B  C  V(B°)  (see  Avellar  and  Hale  [1]). 

We  need  the  following  fundamental  result  from  [l]. 


5.1.  Fix  BCR  .  The  following  statements  are  equivalent: 

0  +  u 

(i)  There  is  an  r  C(R  )  with  rationally  independent  components  such 
that  Eq.  (5.4)  is  hyperbolic  at  (r°,B). 


(ii)  There  is  an  r  €0R  )  such  that  Eq.  (5.4) 


r°,B). 

iii)  Eq.  (S.4J  is  hyperbolic  globally  at  B. 

(iv)  1  *  I  (r.B)  dSf  {|y(6)|:detjyl  -  ^  ^e^  J  =  0,9  CRM 


(v)  li 


#  T  rN  _PYk'r  M 

E(p,0,r,B)  -  detll  -  J  B^  k  e  k  I,  0  €R  » 


p€  R 


aid  T(r,B)  -  (p:  3  0  with  E(p,0,r,B)  -  0>,  then  0  t  r(r,B). 


The  same  result  holds  with  hyperbolic  replaced  by  asymptotically  stable 


and  (iv) ,(vj  are  replaced  b 


(Iv)’  I(r,B)  c  [0,1) 

(v)’  r(r,»)  n  (0,»)  •  0. 


4 
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With  g(X,r,A,8)  as  in  Relation  (5.2),  define 

a”(r,A,B)  ■  sup  {ReX^O  :  g(X,r,A,B)  -  0) 

(5.6) 

a+(r,A,B)  *  inf  {ReX^O  :  g(X,r,A,B)  *  0} 

and  define  a  (r,A,B)  =  -  °°,  a+(r,A,B)  =  +  »  if  the  corresponding  set  is 
empty. 

We  need  the  following  result  which  is  stated  without  proof.  The  proof  uses 
some  special  properties  of  characteristic  functions  which  we  have  not  used  before. 
The  reader  can  supply  the  details  following  ideas  from  [4,  Ch.  12]. 

Lemma  5.2.  If  Eq.  (5.1)  is  hyperbolic  (asymptotically  stable)  at  (r,A,B),  then 
the  difference  equation  (5.4)  is  hyperbolic  (asymptotically  stable)  at  (r,B). 

If  Eq.  (5.4)  is  hyperbolic  at  (r,A,B),  then  a"(ar,A,B),  a+(ar,A,B)  are  con¬ 
tinuous  in  a  for  a  >  0. 

One  can  now  generalize  the  results  of  Section  2  to  the  neutral  equation 
(5.1)  (see  also  Datko  [3].) 

Theorem  5.3.  (A,B)  €  Hr  (or  sp  if  and  only  if 

[N  i  N  l 
(I  -  l  B.)”1  I  A.  U  0  (or  <  0) 

k-1  K  k-0  KJ 

(ii)  g(iy,  ar,A,B)r<  o  for  all  y  €  R,  y  j*  0,  a  >  0. 


Proof:  The  proof  is  the  same  as  the  proof  of  Theorem  2.1,  making  use  of 
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Theorem  5.4.  If 


<s-7>  p<x-=1 . *.*.»>  *  d*t[>-(,-CAs>u  "s»'k*)I(*»^  1Vir“  ^“)] 

then  (A,B)  €  H  if  and  only  if 


(Hj)  Eq.  (5.4)  is  hyperbolic  globally  at  B 
(H2)  P(0,1,...,1,A,B)  *  0 
(H3)  P(iy,s1,...,sM,A,B)  +  0 

for_aU  y  €  *,  y  /  0,  |  s .  |  .  1,  j  .  1,2,  ....  M. 

[N  u  l 

(I  -  l  B  )-1J  A.  <  0 

k»l  k  k=0  kJ 

Proof:  Suppose  (A,B)  €  H.  Then  Lemma  5.2  implies  (Hj) .  Lemma  5.1,  part  (iv) 
implies  there  are  constants  5  >  0,  n  >  0  such  that  |det(I-£N  Bvs^kl SYkM) I 

k=l  *  1  M 

—  ^  ^or  llsjl"^l  <  H»  j  *  1,2,  ...,  M.  Therefore,  P  in  Relation  (5.7) 
is  well  defined  for  ||Sj|-l|  <  n.  Theorem  5.3  impUei  (H2) .  If  (H3)  is  not 
satisfied  for  some  {y,sy . . . ,sM) ,  y  *  0,  |sk|  *  1,  choose  9k  6  K  so  that 

Sk  *  exp^®k)*  ~y  \  *  °»  k  =  1,2,... ,M.  If  r^  «=  -8^/y,  then  g(iy,r,A,B)  *  0, 
which  contradicts  the  fact  that  (A,B)  €  H.  Thus,  (A,B)  €  H  implies  (H.)-(H,). 

Conversely,  suppose  (Hp-(H3)  and  (A,B)  %  H.  Then  there  is  an  o  >  0, 

y  €K,  y  +  0,  r°  €  <&*)M  such  that  g(iy,  ar°,A,B)  «  0.  This  contradicts  (Hj) 

and  proves  the  first  part  of  the  theorem. 

N  ,  N 

The  condition  ReX  [(I-J  B.)  \  A.)  <  0  is  equivalent  to  saying  that 

k-1  *  k-0  * 

Eq.  (5.1)  is  asymptotically  stable  at  r  ■  0.  Thus,  the  last  statement  in  the 


theorem  is  true. 
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Our  next  objective  is  to  obtain  an  alternative  characterization  of  H,S 
in  terns  of  ReX(Ag)  as  in  Section  2. 

Lema  5.5.  If  Eq.  (5.4)  is  hyperbolic  globally  at  B,  then,  for  any 
sj,...,sjj,  |S®|  ■  1,  j  -  1,2,. ..,M,  Hypothesis  (Hj)  implies  either 

P(iy,as°,...,as°)  +  0  for  y  €  ®,  |a|  <  | 

or 

P(iy,as®,...,as®)  /  0  for  y  €  F,  |a|  >  | 

If  Eq.  (5.4)  is  Asymptotically  stable  globally  at  B,  only  the  first  alternative 
holds. 

Proof:  From  the  hypotheses  on  Eq.  (5.4),  there  is  a  >  0  such  that  for 

any  0  <  6  <  6  j ,  there  is  an  q  >  0  such  that 

N  Ykl  YkM, 

|det^I-J  ^k^l  •  •  *SM  jl  i  ^ 

for  1-  6  |Sj )  <_  1  ♦  6  ,  j  *  1,2,  ...  ,  M.  Fix  Sj,  Is^l  *1. 

j  *  1,2,  ...  ,  M  and  define  Q(X,a)  =  P(X,a*j,  ...  ,  as|J)  for  a  €  X  €  $. 
Designate  the  zeros  of  this  function  by  a(X) .  Then  a(X)  is  neronorphic, 
defined  on  some  Reimann  surface  over  the  X-plane.  As  |X|  -*■  00 ,  one  must  have 

(N  Ykl 

I  -  l  [<x(X)s® J  ...  [a(X)s°]  1  0. 

Thus,  for  large  |X|, 

o(X)  €  {  p  €  |p|  <_  1  -  6j  }  U  f  p  €  |p|  >1  ♦  *  jl 


Consider  the  curve  in  the  a-plane  defined  by  a(iy) ,  y  €  F. 
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F  or  y  very  large  and  positive,  suppose  a(iy)  €  {  p  €  $:  jp|  >.  1  ♦ 

Then  Hypothesis  (H^)  implies  that  |a(iy)|  >  1  if  y  +  0.  Thus,  |a(iy)|  ^1 

for  all  y.  Consequently,  Q(iy,a)  **  P(iy,aSj,  ...  ,  as^)  +  0  for  |a|<l, 
y  €  F.  If  ct(iy)  €  {  p  €  $:  |p|  <_  1  -  }  for  large  positive  y,  then 

Hypothesis  (Hj)  implies  |a(iy)|  <  1  for  all  y  #  0.  Thus,  |o(iy)j  <_  1 
for  all  y  and  Q(iy,ots®)  }  0  for  all  |a|  >  1.  This  proves  the  lemma. 
Corollary  5.6.  If  (A,B)  €  H,  then 

N  .  N 

ReX[(I  -  l  B.)'1  I  A.}  +  0. 
k=l  *  k=0 

N  .  N 

If  (A,B)  €  S,  then  ReX ( A_)  <  0,  ReX[(I  -  £  B.)_i  l  \]  <  0. 

—  -  u  ^*1  K  k*0 

Proof:  Let  Q(X,s)  ■  P(X,  s,  ...  ,  s).  For  real  s.  Lemma  5.5  implies 
either 


Q(iy,s)  i  0  for  0  <  s  <  1,  y  £F 
or 


Q(iy,s)  +  0  for  1  <  s,  y  €F 
Hypothesis  (H2)  implies  Q(iy,l)  j*  0  for  all  y  £  F. 

Thus,  Q(*,s)  has  no  noots  on  the  imaginary  axis  either  for  0  5.  s  S  1  or 
for  1  <  s.  Since 


N  N 

Q(X,1)  -  det(XI  -  (I  -  l  B.)"1  I  AJ, 

k*l  k«0 


N 


it  follows  that  ReX  [(I  -  Y  B.)“*  J  A.  ]  +  o.  This  proves  the  first  part  of 

k»l  *  k«0 

the  corollary. 

The  last  part  follows  from  the  last  statement  in  Lemma  5.5  and  the  proof 


N 


of  Corollary  2.7. 
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R^rk  5.7.  One  can  have  Eq.  (5.4)  asymptotically  stable  globally  at  B 
and  have  (A,B)  €  H,  (A,B)  €  S.  In  fact,  consider  the  equation 


*5*8)  -g£-[x(t)  ”  c*(t-r)]  «  ax(t) 

with  c  real,  0<c<l,  a>0.  The  characteristic  function  is 

X(l-ce'Xr)  -  a 


which  is  *  0  for  X  »  iy,  y  £  R.  Also,  there  is  a  real  X  >  0  such  that 
X(l-cexp-Xr)  *  a. 

Let  us  now  consider  in  more  detail  the  scalar  n  -order  neutral  equation 
N 


(5.9)  ytn)(t)-  1  b  y(n)(t-<^)  -  l  a  y<n‘»(t)  -  I  l 

k«l  K  j=l  J  k-1  j-1  J 


where  ^  >  0,  bk€  R,  aik  g  *  for  a11  Let  b  *  (bi»  •**  »  bN3  * 


1* 


a  ■  (a..,  j  ■  1,  ...  ,  n,  k  ■  0,  1,  ...  ,  N) ,  w  *  (<»>j,  ...  »^) 
J  " 

this  equation  as  a  system  of  first  order  equations, 

N  N 

4-[x(t)  - 1  B^t-w^i  ■  v(t)  * 

for  x-  (y,y(1),  ....  y(n’13).  we  have 


If  we  write 


Bk  -  (bgJ),  b™  •  bk.  bg5  -  0  otherwise 


0  1  0  ...  0 

0  0  1  ...  0 


0  ...  1 


no 


an-M>  *h-2p  *10 


0 


0 


I. 


I 

0 


ank  n-l,k-*  alk 


,  k  ■  1,2,  ...  ,  N 


The  characteristic  function  is 


(5.10) 


N  -Xu^ 

g(A,u,a,b)  »p0(X,a)  -  \  p^(X,a,b)e 


k*l 

n 


P0(X,a)  =  detfA-Ag]  =  X  ^ajQ' 

Pk(X,a,b)  =  bkXn  +  a^k  X1 


n-j 


The  corresponding  difference  equation  is 


N 


(5.11) 

with  characteristic  function 

(5.12) 


y(t)  -  l  b. y(t-u.)  •  0 
k=l  K  K 


N  -Xu. 

e(X,u,b)  *  1  -  l  b.e 


k»l 


One  can  now  prove  the  following  generalization  of  the  result  of 
Ziaotovskii  [7]. 

Theorem  5.8.  For  Eq.  (5.9),  (a,b)  €  H  if  and  only  if 


(5.13)  Eq.  (5.11)  is  hyperbolic 

N 

(5.14)  I  a.  +  0 

L.A  ^ 


k-0 


N 


(5. IS)  |p0(iy.*)|  >  ilpkCiy»'a#b) |  for  all  y  €»,  y  i  0. 

The  vector  (a,b)  €  S  if  and  only  if  (5. 14), (5. 15),  and 
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N 

(5.16)  £  IbJ  <  1,  ReX  <  0  if  P0(X.a)  -  0. 

Proof:  From  Theorem  5.4,  condition  (5.14)  is  equivalent  to  Hypothesis  (Hj) • 
Hypothesis  (H^)  is  equivalent  to 

N 

p0(iy,a,b)  ♦  £  ^(iy.a.b)  ^  i1  0 

for  all  y  €  *,  y  +  0,  |sk|  -  1.  k  «  1,2,  ...  ,  N.  But  this  clearly  is 

equivalent  to  (5.15).  Thus,  the  first  part  of  the  theorem  is  true. 

If  (a,b)  €  S,  then  Lemma  5.2  implies  Eq.  5.11  is  asymptotically  stable 

N 

at  b.  Part  (v)  of  Lemma  5.1  implies  this  is  equivalent  to  £  |b,  |  <  1. 

k*l  K 

Corollary  5.6  implies  ReX  <0  if  pQ(X,a)  ■  0.  This  proves  the  theorem. 

We  can  generalize  Theorem  5.8  to  the  case  where  the  delays  are  dependent, 
but  the  result  cannot  be  stated  in  such  a  simple  fashion.  The  proof  is  the 
same  as  before. 

♦  N 

Theorem  5.9.  Consider  again  Eq.  (5.9)  with  *  Yk‘r*  r  €  )  » 

Yk  *  0,  Yk  ■  (Ykl»  •••  .  YyP.  Ykj  10  integer.  With  P0(X,a),  pk(X,a,b) 

defined  as  in  (5.10),  the  vector  (a,b)  €  H  if  I  and  only  if  (5.13) (5.14)  in 
Theorem  5.8  are  satisfied  and 

N  Yk,  Yiju 

(5.17)  pQ(iy,a,b)  ♦  £  pk(iy,a,b)Sj  •••»>,  *° 

k"l 

for  all  y  €F,  y  t  0,  1  sk|  ■  1,  k  ■  1,2,  ...  ,  M. 

The  vector  (a,b)  €  S  if  and  only  if  (5.14) ,(5.17)  and  (5.16)  are  satisfied. 

6.  Scalar  neutral  equations.  For  first  order  scalar  equations 

N  N 

(61)  4^x(t)  ■  \  1  bkx(t'Yk,r)]  •  v(t)  *  \ml  akx(t"Yk,r) 
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where  €R,  k  «  1,2,  ...  ,  N,  the  characterizatir i  of  the  hyperbolic 

and  asymptotically  stable  cones  can  be  specified  in  terns  of  the  properties  of 
the  solutions  of  the  difference  equations 


(6.2) 


(6.3) 


y(t)  -  l  bky(t  -  Yk*r)  •  0 
k-1 

a0z(t)  -  f  «k2(t  -  Yk-r)  •  0 


It  is  the  purpose  of  this  section  to  obtain  such  a  characterization  which 
generalizes  the  results  of  Section  3. 

As  in  Section  3,  let  b  *  (b^,  ...  ,  b^) ,  a  ■  (aQ,a^,  ...  ,  a^) . 


a(8,a)  -  aft  ♦  J  ak  cos  Yk*0 

u  k-i  * 


(6.4) 


0(0, a)  *  l  H  sin  V® 

k*l 

Exactly  as  in  the  proof  of  Theoren  3.1,  one  uses  Theorea  5.9  to  obtain 
the  following  result. 

Theoren  6.1  .  For  Eq.  (6.1) ,  (a,b)  €  H  if  and  only  if 


(6.5) 


Eq.  (6.2)  is  hyperbolic  globally  at  b. 


(6.6) 


(6.7) 


l  \  *  0 

k«0 

Either  o(0,a)  j*  0  or  a(6,a)  »  0,  0(0, a)  «  0. 


and  then  necessarily  aQ  0.  The  vector  (a,b)  €  S  if  and  only  if 


(6.«) 


J,1^ «  *•  *»  ‘ 


and  (6. 6), (6. 7)  are  satisfied. 
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For  the  equation 

N  N 

(6.8)  .  J  jbjXft-r^)]  -  aQx(t)  ♦  ^ajX(t-rk) 


the  above  theorem  has  a  very  sinple  interpretation  for  the  case  of  stability 
globally  in  the  delays.  In  fact,  Eq.  (6.9)  is  stable  globally  in  the  delays 
if  and  only  if 


N 

l  Ibkl  <  1. 

k*l  * 


N 

I  a, 

k-0  K 


<  0, 


\  .W>  i  laol  • 

k*l 


(6.  10) 
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